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Within the framework of the Color Glass Condensate model, we evaluate quark and gluon Generalized Parton 
Distributions (GPDs) and the cross section of Deeply Virtual Compton Scattering (DVCS) in the small-ss 
region. We demonstrate that the DVCS cross section becomes independent of energy in the limit of very small 
Xb, which clearly indicates saturation of the DVCS cross section. Our predictions for the GPDs and the DVCS 
cross section at high-energies can be tested at the future Electron-Ion Collider and in ultra-peripheral nucleus- 
nucleus collisions at the LHC. 
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I. INTRODUCTION 



During the last decade hard exclusive reactions, such as Deeply Virtual Compton Scattering ( DV CS) , 1*(q) + P — » iW) + p'. 
have been a subject of intensive theoretical and experimental studies 0lBSSlll3.ll.IEI3.Q3IiilIi2.Ill.Il3.Il3l- A particular 
interest has been attached to the generalized Bjorken kinematics, 
Oh' 

> 

(N 
^1" 

O 
oo 



-q 2 


= Q 2 lai 


ge, 




w 2 


= (P + q) 2 


large 


) 




Q 2 


const , 




xb 


2P-q 




t = 


A 2 = (P' - 


P) 2 < 


Q 2 



(1) 



where q is the momentum of the virtual photon; P is the initial momentum of the target hadron; P' is the final momentum of the 
target, and t is the momentum transfer. 

In this kinematics the DVCS amplitude is factorized 0.IE1 into the convolution of the perturbative coefficient function with 
nonperturbative Generalized Parton Distributions (GPDs) of the target. Recently, the leading-twist dominance (validity of the 
collinear QCD factorization) in DVCS on the proton target was demonstrated by the Hall A collaboration at Jefferson Labora- 
tory 03], already at rather low values of Q 2 , 1.5 < Q 2 < 2.3 GeV 2 . 
k> , However it turns out that in experiments with nuclei the virtuality Q 2 is not always very large, and one cannot say how accurate 

j^ ■ the predictions based on factorization are, or, in other words, how large the higher-twist corrections are. One of the examples, 
where this approach cannot be applied, is DVCS on the nuclei measured by HERMES collaboration in DESY IU7I1 . Due to 
small-Q 2 ~ 1 — 2 GeV 2 one has to use other effective models, e.g. Generalized Vector Meson Dominance model (GVMD) IU8I1 . 

At very large W 2 (very small values of Xb), the perturbative collinear factorization is expected to break down due to high 
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densities of the partons Ill9ll . Even for relatively large values of Q 2 when the running coupling constant of the strong interactions 
a s (Q 2 ) is small, the effective expansion parameter a s (Q 2 )g(x, Q 2 ), where g(x, Q 2 ) is the gluon density in the target, becomes 
large. This invalidates the perturbative expansion leading to the collinear factorization. Since in heavy nuclei the parton densities 
are enhanced by the atomic number A compared to those in the nucleon, the onset of the effects associated with high parton 
densities may take place at the values of xb, which will be already achieved at the future Electron-Ion Collider (EIC). 

In this paper we use the framework of the Color Glass Condensate (CGC) model offered in J20l 12111 (see also recent re- 
views J22ll23ll24ll25ll ). We generalize the formalism of the CGC model to exclusive reactions and evaluate Generalized Parton 
Distributions (GPDs) and the DVCS amplitude at small-xs. We find that for DVCS off heavy nuclei, the DVCS cross section is 
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virtually xb -independent, i.e. the DVCS cross section saturates in the small-xs limit. The general saturation property built-in 
into this model is an essentially nonperturbative effect, which complies with the general Froissart (unitarity) bound fl26l l27l 12811 

F 2 (x,Q 2 )<ln n (s)^W l (^\, (2) 

where n — 3 for DIS on nucleons; n = 1 for DIS on heavy nuclear targets. This should be compared to the Froissart bound for 
the case of hadron-hadron scattering, a < In 2 s. 

For comparison, estimates of the x-dependence based on the perturbative evolution equations do not possess saturation: The 
DGLAP predicts a fast growing x-dependence 11291 13011 
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xg(x, Q 2 ) ~exp [ 4^ / -a s (Q 2 ) In ^ In- , (3) 



and the BFKL framework 122113 UI32LI33L 13411 predicts the power-growing x-dependence 

xg(x,Q 2 )~[-\ . (4) 

The crucial parameter of the CGC model is the saturation scale Q 2 (x, A), which gives the threshold for transition to saturation 
regime. The saturation scale Q 2 S (x, A) comes into play as a universal parameter in many tasks. For example, in CGC explanation 
of the geometric scaling |3al in DIS data from HERA, the structure function F%{x, Q 2 ) is represented as a function of only one 

variable, i.e. F 2 (x, Q 2 ) = f 



, JIM, 
The paper is organized as follows. In Section ITTA"! we give a brief overview of the model used for evaluations. In particular, 

we generalize the original framework of 120.12111 to the finite nucleus case in order to consider the off-forward matrix elements. 

In Section [HI] we evaluate the quark GPDs, and in Section [TV] we evaluate the DVCS amplitude. In Section \V\ we present our 

results and draw conclusions. 

II. GENERALIZED PARTON DISTRIBUTIONS IN THE COLOR GLASS CONDENSATE MODEL 
A. Overview of the Color Glass Condensate model 

The basic assumption of CGC is that one can separate the partons into fast (xb ~ 1) and slow (xb <C 1) ones, according 
to their light-cone fraction p + . The former are considered as classical "sources", and the latter are the dynamical degrees of 
freedom in the model. In the leading order over a s (Q 2 ) one has just ordinary Yang-Mills equations for the gluon fields, in NLO 
one has a standard loop expansion. It is assumed that dynamics of the "fast" partons does not depend on the "slow" partons; 
thus the configurations of the fast partons are random and one must average over all possible configurations of these "sources" 
Ju(x) = 5^+p a (x), where a is a color index, and x is a coordinate. The weight functional W[p] encodes the dynamics of the 
"fast" subsystem and comes as an external parameter in the model. There are no restrictions on this functional except for the 
obvious gauge and Lorentz invariance. An additional requirement of color neutrality, 



d d x^ a (aV (0)^=0, (5) 

was introduced in Il36ll . It reflects the fact that the physical states are colorless. 

If we define xa as a scale which separates "fast" and "slow" partons, then the dependence of the functional W[p] on the scale 
xq will be described by a kind of "renormgroup equation" 

' d^yiriT^x & v) ttt^ w Ip; t 1 . (6) 



8t 2 J a Sp a (x) /XK ,a ' Sp b (y) 

where r = In ( ^- 1 and % (S, y) is a complicated functional of the field p. 

While in the general case this equation has not been solved so far, there are known solutions for some special (asymptotic) 
cases. Conventionally W[p] is chosen in a Gaussian form J20l |2U 12311 
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where J\f is the normalization factor fixed from the condition J T>pW[p] = 1 and the function X(x, y) is either a constant or 
a function fixed with some additional assumptions. Physically the function \(x, y) describes correlation of partons inside the 
target. It is obvious that in the infinite nuclear matter it may depend only on the relative distance, i.e. 

\{x,y) = n 2 A {x-y). (8) 

In the general case, the shape of the function fJb\{r) is unknown. However, the color neutrality condition (0 and the re- 
quirement that in the low parton density limit the model should reproduce BFKL predictions ©, fix the short-distance and 
large-distance behaviour. It was proposed in [230 that one can use the interpolation 
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Parameterization I: „i(jri = / ^'dK'" = J ^ e V TTtM^ ' (9) 
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where 7 = L/l + jc(3) ~ 0-644 is a numerical coefficient. 

There are also simpler versions of the model |20ll2lll . which neglect correlation of the partons, i.e. 

1 S(x — y) 

Parameterization II: = — — , (10) 

\{x,y) A A (x ) 

where Xa(x~) is some function 047I1 . In subsequent sections we will consider first evaluation with a simple parameteriza- 
tion ( [Tol l, and after that discuss, how the results change for the parameterization (|9). 

The choice of the Gaussian parameterization ^} enables us to evaluate all the results analytically. Notice however, that (j7]) is 
explicitly C-even, i.e. the number of quarks is equal to the number of antiquarks inside any target in this model. This agrees 
with experimental fact that the quark and anti-quark parton densities are approximately equal at small xb- On the other hand, 
C-parity of (0 implies that the model does not distinguish matter and antimatter and is not applicable to evaluation of some 
quantities. For example, the baryon number and electric charge of the target are exactly zero, since they are due to the valence 
quarks. 

An interesting generalization of the Gaussian parameterization Q was discussed in H38f1 - In particular, it was found that 
for the model of k 3> 1 independent noninteracting quarks the distribution is indeed Gaussian, and the first correction is 
proportional to ~ d a bc f d 3 x p a (x)p b (x)p c (x), where d a bc is defined from the anticommutator of the generators T a of the 
group, {T a , T b } = 2d a b c T c . However, for the DVCS and singlet GPDs discussed in this paper the C-odd correction does not 
contribute. 

It is well-known that at high-energies, the real par t of scattering amplitudes is suppressed by the slow energy dependence 
of the amplitude compared to the imaginary part 11391 l40l 14111 . Therefore, it is sufficient to consider only the imaginary part. 
Actually, as we shall show in Sect. IIV1 the real part of the DVCS amplitude in the CGC model is exactly zero. 

The generating functional of the model has a form (480 

zm = jD P w W s { DM { A+)e ^ P] . (id 

where S[A, p] = S[A] + J dxp a (x)A c !_{x) and we used light-cone gauge n ■ A = 0, n 2 = 0. In order to restore the explicit 
gauge invariance of the action S[A, p], the interaction term J dxp a (x)A^ (x) is sometimes replaced with Tr J d 3 xp(x)W[A, x], 
where 

W[A, x\ = Pexp ( ig f d(A + {() J (12) 

is the Wilson link. 

B. Finite nucleus 

Since in this paper we are interested in DVCS-off-forward reaction, we can no longer use the infinite nuclear matter approx- 
imation. Indeed, the DVCS cross section off a nuclear target rapidly decreases as one increases the momentum transfer t. As a 



result, the sizable cross-sections exist only for |i| ~ \jR\, where R A is the nuclear radius. In the infinite nuclear matter, all the 
off- forward cross-sections vanish 14911 . This means that we have to take into account the off- forward kinematics from the very 
beginning. If the coordinate of the nucleus center of mass is X, then the weight functional W[p] may be chosen as 

1 fja„anx v , , d ,P a (x-X)p a (x-Xy 



W p [ Pl X] = c* v l--J d*x6(\3! ± - X ± \ < R A Y \,__ K x . > , d3) 



where we extracted the "zero mode" (integration over the nucleus center of mass) explicitly according to standard technique 142 
and introduced an explicit cutoff factor 6(\x± — X±\ < Ra) which forbids the color condensate p a {x) from outside of the 
nucleus. The cutoff in x~ is provided by the factor \ A (x~~ — X~). The interaction of gluons with the condensate is also 
modified by this cutoff factor: 

S[A,p] = S[A] + Tr f d 3 x9(\x ± - X ± \ < R A )p(x)A_(x) (14) 

for the linear interaction, or 

S[A, P ] = S[A]+Tr f d 3 x6(\x±\ < R A )p(x)W[A](x) (15) 

for the interaction via Wilson link ( TTZl i. The generating functional (fTTT i takes the form 

f ,-Xy r , [DA5(A+)e iS l A >P)-f dx i- A 

Z\J] = J DpdXe^W[p,X] J j K DM [ A+)eiS{A , p] • (16) 

Notice that the formal introduction of the ^-functions is equivalent to the redefinition of the functional integral: 

I Dp := IJdp(a?-, \2± -Xx\< R A )d 3 X. (17) 

Indeed, configurations with /o(|^j_| > R A ) ^ do not interact with anything and thus contribute only to the normalization 
constant. 

Since the coupling constant a s is small, we can take the integral over the gluon field A p in dTTb in the saddle-point approxi- 
mation. In the leading order, the gluon field A^ is just the solution of the equation of motion 

D,F^(x) = S v '+S(x~)p a (x ± ), (18) 

where p a (x±) is the arbitrary external field, and an additional gauge constraint A + = is implied. Notice that we do not 
impose any conditions onto the gluonic fields A a at the large distance |af| > R A . The solution of the equation ( fT8l is 12311 

A" = U (A„ + -d^ U\ (19) 



where 15 

A„ = S fi+ a(x-,x±), (20) 

U = Pcxplig dz-a a (z-,xj_)T a I , (21) 

1 f 1 1 
a(x~,x±) = — xrP= d 2 y ± — In- —2 p{x~ ,y ± ) , (22) 

p(x~,xx) = U J< (x^ 1 x±)p{x^ 1 x 1 _)U(x^,x±), (23) 

and T a are the generators of the color group. 

The analytical solution < fl~9T > enables us to evaluate different correlators. A straightforward evaluation of the (pp)-correlator 
with the weight function ( fT3l yields 



P' \p(x)p(y)\ P)=P + d 3 Xe tAX 9(\x ± - X ± \ < R A )X A ( X - - X-)5 3 (x-y) = f(A)P+e- lAx S 3 (x-y) , (24) 
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We can see that for any fixed nonzero Aj_ 7^ the result vanishes in the Ra — * 00 limit in agreement with discussion at the 
beginning of this section. The evaluation of the gluonic GPDs defined as [fllTl 
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(27) 



is done in quasiclassical approximation, 

xH3(x,£,t)* Jd 3 Xe^ J dz-e^ 1 '- F$ k I -^- - X \ F? a I ~— - X 1 , (28) 

where Fn V in the rhs of j28] i corresponds to the classical solution found in previous subsection. Evaluation of ( f28b [23] gives Iplll 
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and j a (x , f± ) is defined as 
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(29) 



(30) 



d 3 k e -ix(k+A/2) iy(k-A/2) 



(31) 



As one can see from ( |29l , the gluon GPD H 9 (x, £, t) has a trivial x-dependence 1/x for all (£, £), since x does not enter the 
right-hand side of Eq. ( 1291 . Physically, the exponent in d29l takes into account nonlinear in a s effects in the model. 



C. Alternative kernel 



In this section we discuss how all the previous formulae change with an alternative weight function (|9). The weight functional 
in this case should be written as 



W p [p, X]=Afcxp! ~ I 



x I - X J 
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(32) 



where the function p A (z) describes correlation of hadrons inside the nuclei and was defined in ([9). Performing the evaluation as 
was discussed in Sect. Ill Bl we obtain 
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(36) 



From d33l l we can see that in the finite nuclei the color neutrality condition (01 implies that we have to identify p? A from (O 
with p 2 A {r). For ja {x~ , f±) we can immediately obtain 



(P'\a(x)a(y)\P) 



d 3 k 

WT 3 



£i(*)- 



e -ia:(fc+A/2) e iy(fc-A/2) 



*(*- - y-)7A(x-,i?± - fae^W* . (37) 



Thus we can see that this kernel differs from the previous one only by an additional factor p 2 (k) in the integrand. 

D. Quark propagator in CGC field 

Although for the evaluation of the DVCS amplitude one may use the color dipole approximation, in this paper we evaluate the 
GPDs and Compton amplitudes directly. In the diagrammatic language this corresponds to summation of all J52I1 the multigluon 
diagrams, whereas the color dipole approach assumes either only Born term contribution or Eikonal approximation, as is shown 
on the Fig. Q] 




Figure 1: Diagrams contributing to DVCS in color dipole approximation. See 14311 for an example of DIS evaluation in this approach. 

For the evaluation of the quark GPDs in the leading order over a s (Q 2 ), we need to evaluate the quark propagator in the 
classical gluonic field found in the previous section. To this end, we consider only the zero width limit, 



P = 5(x )p(xj_). 



(38) 



Beyond this limit, equations with explicit x~ -dependence become much more complicated. Physically, the use of d38l l in the 
off-forward kinematics is justified, since the light-cone fractions of the partons are small, i.e. i,^<1, 

The basic idea is that for x~ ^ the field p(x) = and we have just vacuum equations, gluon field A^ reduces to a pure 
gauge. It is possible to choose the gauge in such a way that for x~ < the field disappears, A^ — 0, and for x~ > it is a pure 
gauge, A^ = -UdJJ* and thus the wave function of the quark has a form 



~>Pps{x) 



u s (p)e~ ipx , x~ < 

Jd 4 p'5(p 2 -p' 2 )j: s , C ss ,(p,p')u s ,(p')e- wx ,^ > 



(39) 



where u s (p) is a free Dirac spinor, and the matrix C SS ' (p, p') is found from the continuity at the point x =0. One subtle point 
is that the Dirac operator has the form 



iD = id- j 



(40) 



and matrix 7 is singular, because it is proportional to the light-cone projector A' ', This implies that the continuity condition 
must be imposed not on the function tp ps {x) as a whole, as in J20ll2ll 12311 . but rather only on the component 15311 ip ps (x) = 
^ ^Pps(x). The final result for the wave function is 

^ ps {x) = 6{-x-)u s (p)e- ipx 

+6(x-)U(x ± ) f^Ls(k--p-)S (k+ ^ +p2 



{2-Kf 



2p- 



x (fd i ze i ^- k ^- z U\z)\e- ikx (l + -^y=0e± + M) ) A<->u»(p), 



where M is the mass of the quark. The evaluation of the quark propagator according to 

d4 p E s ipps(x)i>ps(y) 



S(x,y) = 



(2?r) 4 p 2 - M 2 + tO 



yields 
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jd 2 ze l( - p ^ k ^ z (l + ^^(h + M)) A(-\p + M) (U{x ± )W (z) - 1) , x~ > 0,y~ <0 
where 5o(x — y) is the free propagator 

d 4 p e^P^"^ 



Sbfc- y) = 



(2tt) 4 p - M + tO ' 



(41) 



(42) 



(43) 



(44) 



It might be checked that the propagator S(x,y) satisfies the equation (iD — M)S(x,y) = 5(x — y) as well as reduces to 
Sq(x — y) in the U — ► 1 limit. 



III. UNINTEGRATED QUARK GPDS 

In this section we evaluate unintegrated quark GPDs defined via the following matrix element (we assume that the target has 
spin 0) 



H(x,^Aj_,k A 



(45) 



t / ** 



ikj_rj_ ixP^~ z 



p' 



*|-T-ThMV + T 



In the forward limit (Aj_ — > 0, £ — > 0) the function # I x, £, Aj_ , A:^ ) reduces to unintegrated parton distribution q[x,k±\, 
and when integrated over k±, it gives ordinary GPDs. In the quasiclassical approximation, d45T > reduces to 

H(x,t,& x ,iex)= (46) 



f^L_ e ixP+z- f d 2 r±e -^r ±i p+ f^x^AX / Tr 



■''_ 



r± 



' \ 2 2 '22 

where here and below angular brackets without explicit initial and final states (...) are the short-hand notation for averaging 
(integration) over all possible configurations p(x), i.e. (0):= J VpW[p]0(p). Substituting the propagator d43b and taking 
the integral over each domain, one obtains the final result 



H (x,£,Ax,kj_) = #<+-) + # ( - +) 



(47) 
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the superscript signs (±, ±) refer to different integration domains over (dX ,dz ) in j46l ), and function 7 ( k — ^-, k + 4p 



is defined as 



A x , A 1 



J^ e ^^ p{u ^ x + I )u{x 



Evaluation of this quantity (see Sect|B]for details) yields 






d 2 re^ F d 2 X,e iA ^ 



'f(d,£-X)+f[0,-Z-X] 



(50) 



(51) 



cxp 



Notice that GPD ( f47T > is antisymmetric, i.e. H(—x, £) = —H(—x, £). Also, ( |47] > is not required to satisfy polynomiality since 
the original model is valid only for i<l. 

One of the subtle points of the result (l47b is the logarithmic behaviour ~ In \x ± £| in the vicinity of the points x ~ ±£. 
Physically, in these points one of the quarks has a zero light-cone fraction and becomes especially sensitive to the details of the 
model. However, since we are in a saturation regime, factorization formula does not work and we expect that such behaviour 
should not cause any physical problems. In Appendix lAl we give details of evaluation of d47b . and in particular discuss the 
logarithmic singularities. 



IV. DVCS AMPLITUDE 



In this section we evaluate the DVCS amplitude directly (not using factorization). The first reason for this is that the GPDs 
evaluated in the previous section are valid only for the small x <C 1 whereas the convolution formula which follows from 
factorization implies integration over the light-cone fraction over the region — 1 < x < 1. The second reason is that, as we 
discussed in Sect. ©, in the saturation (high-density) regime the convolution formula becomes invalid. 

The starting point of our derivation is the definition of the DVCS amplitude 



(52) 



In the quasiclassical approximation the matrix element (P' \J v (0)J^(z)\ P) A is reduced to 



(P' \M0)Mz)\ P) A = J d 3 Xe tAX (P' J u (-X)J^z - X) 
rf 3 Ie tAi (Tr [ llx S{z - X, -X)^S(-X, z - X)\) 



P 



(53) 
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where S(x, y) is the propagator d43b . Substituting (T43b into d53l ) and taking the integrals, we may reduce the DVCS amplitude 
to the form 



Afiv — i 
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x 1 8N C S M .S U+ (m 2 -%l- A) Ul 2 - f ± + ^ J 




+8N C S. ± S V± 



Pip! {{q-f - 4(p-) 2 ) +p-g MV (4p~(M 2 +pl) - q-p ± ■ Al) 



(^) 2 -^^ 



^2 (9 ) \ f i/2 r,^2 L2 , A i 



+ 87V c (^+^_ + <W V) I (P") - ^^H I 2^/ 2 - 2pl - *i + -^ I I . (54) 

One interesting point is that the real part of d54b is exactly zero. Indeed, taking the imaginary part of the first ratio containing 
— iO and using 

—^— = p(-)+inS(x), (55) 

x — iv \xj 

we can immediately find that the argument of <5-function is zero only for 



q~ 2(kj + AP) q- 

2 V + Q 2 " 2 



|P-| = V^/1+ V± o 2 L >- (56) 



i.e. outside the integration domain. For comparison, from phenomenology it is known that the high-energy amplitude gets 
dominant contribution from the imaginary part. 

V. RESULTS FOR GPDS AND DVCS CROSS-SECTIONS 

In this section we present results of the numerical evaluation of the GPDs and DVCS cross-sections. In subsection (IV Ai we 
consider first the results with a simpler parameterization ( [Tol l, and after that in subsection flVBI i with a more realistic parameter- 
ization (0. 

A. Results with parameterization II 

As one can see from d29| i, for both parameterizations I and II the x-dependence of the gluon GPD H A (x, £, t) is trivial-just 
1/x for all (£, t). For quark GPDs H A (x, £, t) the x-dependence is more complicated, however in the forward case the parton 
distribution qA{x) has also a simple 1/x-dependence. For better understanding, we prefer to discuss out results for the gluons 
in terms of the ratio H A (x, £, t)/gA(x), which measures the off-forward effects, and gA(x) is the forward gluon PDF evaluated 
in the same model. 

In Figure|2]we plot the £ and ^-dependence of the ratio H A (x, £, t)/gA(x) for different nuclei. In the left panel of Figure|2] 
we plot the ^-dependence of the ratio H A (x, £, t)/gA(x) in nuclei for £ = 0. We can see that H 9 (x, £, t) is decreasing as a 
function of t. For the sake of comparison, on the same plot we also plotted in grey lines the nuclear form factors in conventional 

exponential parameterization, .FU(t) = cxp I -g 1 *), and for radius Ra we used Ra = 1.2 fm x A 1 / 3 . We can see that to a 
good extent the i-dependence of the GPDs is similar to that of the form factors. 

In the right panel of Figure [2] we plot the ^-dependence of the gluon GPDs in nuclei. We can see that in the small-£ region 
H 9 (x,£,t) is independent of the skewedness £, This results is quite easy to understand: in the ultrarelativistic limit the nucleus 
in laboratory frame is squeezed to an infinitely thin "pancake", so the condensate distribution along the :r~-axis is strongly 
peaked around x~ w 0, Xa(x~) ~ S(x~). As a consequence, the gluon GPD which is proportional to the Fourier of Aa(x~), 
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almost does not depend on A + ~ £, The only exception is the region of sufficiently large £ ~ 0.1, where the ^-dependence is 
mainly a kinematical effect-the increase of H 9 (x,t;, t) is due to decreasing Aj_ at fixed t. However, these values of £ ~ 0.1 are 
too large, and our extrapolation of the model becomes unreliable. 



H' A (x,£,t)/g A (x), Q 2 = l GeV 2 ,%=0 



W (x&t)/g (x), Q-=I GeV, t=-0.01 GeV' 
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Figure 2: Left plot: ^-dependence of the gluon distribution for different nuclei. £ = 0, Q =1 GeV . For comparison, we also plotted in grey 

lines the nuclear formfactor in the simplest exponential parameterization Fa it) = exp ( -£-t ) . Right plot: ^-dependence for the same nuclei 

for fixed t = —0.01 GeV 2 , Q 2 = 1 GeV 2 . We do not plot the a- -dependence of the gluon GPD H 9 , which is according to {29} just a trivial 
1/a; for all (£,£). 

In Figure [3]we plot the x-, £- and t-dependence of the quark GPD H A (x, £, t) in nuclei. As one can see from d47b . in the 
forward limit the quark distributions have a very simple ^-dependence, H A {x, 0, 0) = q A (x) ~ 1/x. For better legibility, we 
prefer to discuss out results for the quarks in terms of the ratio Ha (x, £, t)/qA (x) , which measures off-forward effects. 

From the left panel in Figure[3]we can see that for x -C £ the GPD Ha(x, £, t) is decreasing approximately as Ha{x, £, t) ~ x 



and as a result the ratio HA(x,Z,,t) / q A (x) behaves approximately as H A {x ,(,,€) /qA{x) 



For a- > £, H A (x,£,t) 



qA{x)FA{t), and the ratio is a constant. In the point x — £ we have a singularity ~ In |x — £|, which was mentioned at the end 
of Sectionllllland discussed in details in AppendixlAl 

From the middle panel in Figure[3]we can see that as a function of £ the generalized quark distribution is a constant for £ <C x, 
but is a decreasing function for £ 3> x. 

From the right panel in Figure |3]we can see the t-dependence of the GPD Ha(x, £, t). For the sake of comparison, on the 
same plot we also plotted in grey lines the nuclear form factors in the frequently used exponential parameterization, Fa (i) = 

exp ( -gM J . We can see that Ha(x, £, t) is decreasing a bit faster than Fa(£). 
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Figure 3: Left plot: ^-dependence of the quark GPD H A (x, £,£). £ = 10 2 , Q 2 — 1 GeV 2 . Middle plot: ^-dependence of the same GPD 
H A (x,£,t) for fixeda; = 10~ 4 ,Q 2 = lGeV 2 £ max = y/-t/(4M 2 - t). Right plot: t-dependence of the same GPD H A {x, £,£). For 
comparison, we also plotted in grey lines the nuclear formfactor in the the simplest exponential parameterization F A (t) = exp ( -^-t ) . 

In Figure|4]we plot the £- and t-dependence of the differential DVCS cross-section da/dt for fixed Q 2 and different nuclei. 
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From the left part of Figure |4] we can see that the cross-section is growing when £ is decreasing, but at some £, which we call 
£,sat{Q 2 : A), we have a qualitative transition to the saturation. The value t; sa t{Q 2 , A) depends on the external kinematics. The 
relatively large value £ sat ~ 0.01 is due to the small value of Q 2 — 1 GcV , £ SQ t is decreasing when Q 2 increases. 

From the right plot on Figure [4] we can see the i-dependence of the differential cross-section da/dt. For the sake of com- 
parison, on the same plot we also plotted in grey lines the nuclear form factors in conventional exponential parameterization, 

6 



Fa{1) — CX P ( -ftt) ■ We can see that da/dt is decreasing a bit faster than F\(t). 
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Figure 4: Left plot: ^-dependence of the differential DVCS cross-section in the CGC model for different nuclei. Kinematic is chosen as 
Q 2 - 1 GeV 2 , t - -0.01 GeV 2 . Right plot: t-dependence of the DVCS cross-section at fixed £ = 10~ 4 . On the right plot, we also plotted 
in grey lines what one would have with the simplest factorized i-dependence of the DVCS amplitude and exponential parameterization for the 

formfactor: ^ ~ F\{i) ~ exp (^t 



B. Results with parameterization I 

In this section we discuss the results of Color Glass Condensate model in parameterization (0. The crucial point is that this 
model explicitly contains the saturation scale Q 2 and as a consequence we can apply it only to the kinematics where saturation 
is present. In our evaluations we used for Q 2 the parameterization from 12411 . where Q 2 (A) is found as a solution of the equation 



](A) = a s (Q 2 )N c ^ 2 A ln 



Q 2 M) 

. A QCD , 



(57) 



This equation has real solutions only for A > A min (Q 2 ) ~ 150 for Q 2 ~ 1 GcV 2 , and A min (Q 2 ) is a growing function of Q 2 . 
InFigures|5]and|6]weplot the £ and i-dependence of the gluon and quark distributions HA(x,t;,t)/qA(x), and in Figure|7]we 
plot the £- and i-dependence of the differential DVCS cross-section da/dt . We can see that qualitatively the behaviour is the 
same as in the previous section, although absolute values differ. 



C. Comparison to DVCS cross section in GVMD model 



In the Figure [8] we compare predictions for the DVCS cross-section with our earlier result II 1811 obtained in Generalized 
Vector Dominance Model (GVMD). We can see the difference in predictions of GVMD and CGC models: In contrast to the 
saturation behavior in CGC, the GVMD cross-section is slowly growing as £~" when £ is decreasing. Nevertheless in the region 
10 -5 < £ < 10~ 3 predictions of both models have comparable values. 



VI. CONCLUSION 



In this paper we considered Generalized Parton Distributions (GPDs) and Deeply Virtual Compton Scattering (DVCS) ampli- 
tudes in the Color Glass Condensate model. We modified the original formulation of J20l 12111 to off-forward kinematics of hard 
exclusive reactions, which provided the necessary framework for the calculation of GPDs and the DVCS amplitude. 
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Figure 5: Left plot: ^-dependence of the gluon distribution for different nuclei. £ = 0, Q — 1 GeV . Right plot: ^-dependence for the same 
nuclei for fixed t — —0.01 GeV 2 , Q 2 — 1 GeV 2 . We do not plot the x-dependence of the gluon GPD H 9 , which is according to ( |29t just a 
trivial 1/x for all (£, t). On the left plot, we also plotted in grey lines the nuclear formfactor in the the simplest exponential parameterization 

J FA(t)=e X p(^i). 
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Figure 6: Left plot: ^-dependence of the quark GPD Ha(x, £, t). £, — 10 , Q — 1 GeV . Middle plot: ^-dependence of the same GPD 
Ha(x, £, t) for fixed x = 10 -3 , Q 2 = 1 GeV 2 . Right plot: t-dependence of the same GPD Ha{x, £, t). On the left plot, we also plotted in 

grey lines the nuclear formfactor in the the simplest exponential parameterization FA(t) = exp ( -g^i J . 



We evaluated the quark and gluon GPDs in this model and studied their dependence on variables x, £, t. We found that the 
gluon GPD H 9 in this model has a simple ^-dependence H 9 (x) ~ 1/x for all (£, t). Similar 1/x-behaviour was observed for 
the quark GPDs H q in the x ^ £ region and in the forward limit (t = 0). Both the quark and gluon GPDs are decreasing as a 
function of momentum transfer t, and the quark GPD is decreasing a bit faster than gluon GPD. 

Without assuming the validity of the collinear factorization, we evaluated the DVCS cross-sections in the small-£ region on 
the large nuclei. We found that in this region the DVCS cross-sections are almost independent of £. This is a manifestation of the 
general saturation property inherent to the CGC model. As far as absolute values are concerned, we found that the predictions of 
CGC in the relevant range of £ are comparable with predictions of other models, e.g. GVMD. Currently there is no experimental 
data available for DVCS cross-section in this kinematics. 

The present calculation should be important for a wide range of the future experiments. For example, gluon GPDs in the small- 
x region may be used for evaluation of the heavy vector meson production in ultraperipheral collisions at the LHC 1451 14611 . 
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Appendix A: DETAILS OF EVALUATION OF 06) 

In this section we evaluate the unintegrated GPD j45] l, which in quasiclassical approximation was reduced to 



H(x,£,A ± ,k ± 
dz~ 



iP + 
iP + 



#r\e- ik±fl 



d 3 Xe- iAX (Tr 



1 + S 



r -±-x — 

2 ' 2 



r± 



X 



(2tt) 3 (2tt) s 



^•fe-ipi-fi (Tr [7+5(6; 6)]), 



(Al) 
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where we changed the integration variables according to 



and introduced shorthand notations 



U 2 2 



t - Z -u Tl - 

6 - T + T 



x. 



A 



pi = a-P + + /cj_ - — . 



p 2 = a;P + + fcx + — 



(A2) 
(A3) 

(A4) 
(A5) 



Now we have to consider separately the first case £j~ > 0, £% < and the second case £-f < 0, ^J" > 0. All the other regions 
are just the vacuum contributions ~ 5 2 (Aj_) and must be omitted. For the sake of brevity we will refer to the contribution of 
the first region as H^ , and to the second one as H ^. 

For ^-functions of arguments ±£12 we will use an integral representation 



o(±0 



2ni 



da- 



a — i0 



1 

2tt£ 



da- 



a + i0 



(A6) 



1. Evaluation of H A 



In the first case we have explicitly 



, d 3 f^ d 3 fn ■ , f dn-idno ^("iSi -"262 ) 



d 4 p 



(2tt) 3 (2tt) 3 
1 /■ A 



(2tt) 4 p 2 - M 2 + i0 J (2tt) 



cxp < —z 



(2tt) 2 (ai-i0)(a 2 -i0) 
gi + M 2 



2p- 



"£l - 3-L ' £l-L - P + £ 2 + P-L^ 



x / d a «e i(p - L -« J - ) *Tr 
Now evaluate each of the integrals: 



7+ (l + -=j=($± + M)) A<-'(p + M) 



(f/t(z)[/(6x)) 



(A7) 



d£id£. 



(2ny 



'i-piPlZl -W2 S2 pif^iCi -"2^2 ) , 



.?I + M 2 _ . .,_ 

ex P I -* — Si + l P ?2 



aUi+pf- 2 _ J(5(a 2 +P 2 f -P + ) , 
Now change the dummy integration variables £/- and z to X 1 - and /r 1 : 






_ 



tf-^-V 



z :=X 



(A8) 



(A9) 



(A10) 



fd 2 ze^-^ z J d ^p e -ipUt+i P i& e i«±St-iP2±& (Ufa±)U*(z)) 



(27T) 



(All) 
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where 7 was defined in ( T50b . It is convenient to make a shift of the dummy integration variable according to 

where k = k — q . 



(2tt) 2 J (2tt) 



(A12) 



=> H + ~(x,£,t,k±) = -i/V, 
1 



2^-L 



d k 



A± JL AJ 



M 2 - p 2 ±(fcl - K±) 



(2tt) 2 ' V" ' 2 ' 2 



dp + dp 



(2tt) 2 2p+p- - p^ 2 - M 2 + iO 



P+-P+- iO (^-«x) 2 +a/ 2 _ p + _ - p - 

Now we take the integrals overp + ,p~ in ( IA13b . The first integral is taken overp + , the result is 



(A13) 



dp* 



1 



1 



i6(p~) 



(2?r) 2p+p~ - p^ 2 - M 2 + iOp+ -p+ - iO 2p+p~ - (p 4 -) 2 - M 2 + iO 



Integration over p yields 

H+-(x,Z,t,kj 



2N r 



d 2 K A 



A J 



(2^ 7 r + ~'^ 



A J 



M 2 -p 2 ±-(k±_ 



(x - £)((fcjL - «±) 2 + ^ 2 ) - (* + 0(P2 2 + M 2 ) 



In 



a; - £ ((*x - *x) 2 + AP) 



+ e ipr + m 2 ) 



2N r 



x A± 
7 I « +-S-.K" 1 " - 



A- 



d 2 K J 



(2tt) 2 ' V" ' 2 ' 2 



M 2 - ( fc + %H • ( A; - k_l 



(x-0 ( (fc-^V + M 2 J -(» + £) U^+%V 



In 



z-£ 


(J- 


Kj_j 


+ 71/ 2 


x + £ 


(fc + 


*)" 


+ M 2 



(A14) 



(A15) 



2. Evaluation of H 



In complete analogy we evaluate the term H + : 



(2tt) 3 (2tt) 3 
d 4 p 1 /• d 2 g 



J (2ir)° {2tt)° 
f d 4 p 1 /• d 2 . 

X 7 (27r)i p 2 - M 2 + iO J J^) 2 ' 

I d 2 ze- i(p± - q±)z Tr 



(2tt) 2 (ai + i0){a 2 + iO) 



2p- 



Now take the integrals term-by-term in complete analogy with the previous case 



/ ^^e*f «r-*f ff e^xC—ff) cxp (i^^e 2 - - ip + (I 

q\ + M' 2 



(UHtii)U(z)) 



5 («i +pf — p + ) 6 I a 2 + p\ 



2 P - 



J (27TJ 2 J J (27T) 4 



(A16) 



(A17) 



(A18) 



16 



Now change the dummy integration variables £^-, z to X, p 1 - according to 



2 ' 



x- 



(A19) 
(A20) 



^ fd 2 ze- l ^-^ z f d2 l^y e -*ttt+tf& e -H±£+ie±tt (uH^)U(z)) 



(2tt) 2 



(2tt)< 



where function 7 was defined in ( l50b . Now shift the dummy integration variable according to 



(2tt) 2 J (2tt) 2 
1 



where k 



d 2 ^ ( , A- 1 - j A 1 

M 2 -P;U_(£l-Kl) 
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dp + dp 



(2tt) 2 2p+p- - pj 12 - M 2 + iO 
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First take the integral over the p + 
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1 
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lV{-p 
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next take the integral over p~ : 
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X (1 + 0((fcj- - «±) 2 + ^ 2 ) - (* - 0((Pi-) 2 + M 2 ) 
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(s+oi (fci-*.ij + A/2 ) - (* - n * 

In summary, we have 
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(A25) 
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Notice that the sum 

H (x, e, t, fcj.) = ff+~ (x, e, t, *x) + H-+ (x, £, t, fc ± ) (A27) 

is antisymmetric w.r.t. the inversion of the light-cone fraction x — > —x, i.e. H I —x, £, t, kx ) = —H (x,£,t,k±) . 

We can see that in the points x = ±£ the result (!A27b has logarithmic divergences ~ In \x =F £|. Physically, in this points one 
of the quarks has a zero light-cone fraction, and as a consequence ( |A27t becomes very sensitive to the details of short-distance 
structure of the model. When we evaluated dA27| ), we integrated over jr 1 up to infinity. Rigorously speaking, this contradicts the 
basic assumptions of the model, in particular, d38l l. which is valid only when the moments of the active partons are much smaller 
than the moment of the whole nucleus. However, since for p^ 2 ^ the integrals were convergent (the dominant contribution 
comes from the region where the model is valid), we could ignore such an explicit cutoffs. Notice that in evaluation of the 
physical DVCS amplitude d54l l the cutoffs \p~ | < q~/2 were provided by the external kinematics. Generalization of (|38l l to the 
more realistic color source is a much more complicated task. 

Appendix B: (U^U) CORRELATOR IN FINITE NUCLEI. 

As we have seen in the previous section, as well as we will see in the next section, physical observables depend on the 
correlator 

(P'\U'<(x)U(y)\P)^P + f d^Xe^Tr^ix-X^iy-X)). (Bl) 



Notice that the weight functional W is expressed in terms of the field p, i.e. the correlator is essentially nonlinear object. 
In the finite nucleus evaluation of this object slightly differs from the original derivation given in luOlulll . However, since the 
weight functional W[p] is Gaussian, the total result can be expressed in terms of the elementary correlator J54ll (P' \pp\ P) . 

The final result of our evaluation is 

S(x,y) = (P'\uHx ± )U(y ± )\P)=e tK ^ 1 ^ f d 2 X e lA ^ x (B2) 

x exp 



';(o,^i-J?)+/(5,-^ 



-x 
~g 2 N c I ^ '—— ± '- -f[2±- y±, -X\ 



where f(r u r 2 ) = J ^e^^ J+~ dz~nf A (z-, n; A). 
Indeed, using definition 



/ f +QC 
U(x) = Pcxp I ig / dz a a (z ,x±)T a 

\ J — oo 



(B3) 



we may notice that 

• Only the even powers of a give nonzero contribution to ( IB2I > 

• The first term (zero order in a) is proportional to 5(A) and vanishes in the off-forward limit. 
Contribution of the second-order term gives 

-g 2 N c e lK ^-^ (cos(a ± X± ~ V± \ Jdz-j A (z-,0 ± )- Jdz-j A (z-,x ± -ij ± ;A)\ . (B4) 

It is very convenient to introduce temporary notation f dz~j A (z~,f±; A) = f(r±; A). In this notation ( IB4b reduces to 

~g 2 N c e^^^ fcos ( ^f=] /(0 ± ; A) - /(r1;A) ) , (B5) 
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where we used notation r = x± — y± . 

Evaluation of the higher-order contributions is a bit more tricky. First we have to notice that the Gaussian form of W[p] 
enables us to introduce a sort of Wick theorem for evaluation of the multileg correlators. After that, we have to make Fourier 
transformation of each correlator, take the integral over d 2 X± and make the Fourier back to coordinate space. Performing such 
procedure step-by-step, contribution of the 2n-th order term after some manipulations may be reduced to 



(-l) n ~ m g 2n N? f d 2 A^ f d 2 A 



2n mm(m,2n— m) 

^ Q ^ k\{m-k)\(2n-m-k)\ J (2n) 2 '"J (2tt) 2 ' 

rm+j«1 



m-j£h^-z« 



IJ f(6,Ai)e lSA 

1=1 

In 

*^E E 

ro=0 fc=0 



*=[E±*]+1 

{-l) n - m g 2n N r c 
k\{m — k)\(2n — m — k)\ 



f(r±,Ai 



vfcl^l+i 



fd 2 Af [d 2 At , ^ KJ 



(27r) 2 '"y (2tt) 2 

n M^y pKt 



n f&&> 



-ir A^/2 



fl f(r±M) 



+1 



Now we replace back 5(&± - £" =0 &t) = J d 2 X e^-^ \[ n l=1 e -^t^± and reduce (|B6) to 



(B6) 



n 



2n Tnin{m,2n— m) 

d 2 Xe lA ^ ^ 

d 2 A^ 



A A- 



n 



_+S± 



(2tt) 2 



f -E E 

m=0 fc=0 

/(0,A>K r ^ /2 
2 f(T±M) 



n—m „2n /\m 



A 7 ,'.' 



fc!(m-fc)!(2n-m-fc)! 

d 2 A 



n 



=[=£=]+! 



(27T) 



/(0,A^)e-- A ^ 2 



(2tt 



2n mm(m,2n-m) 



m=0 fc=0 



(-l) n - m 5 2n Af ( 



A;!(m — fc)!(2n — m — k)\ 



x/ 1 



3, J -*)/[=**]($-£-*) /*(&,-*) 



e iA£ ^ /d 2 X e ^*exp 



'/(o,f-xj+/(o,-§- 



f(r,-x) 



(B7) 



in agreement with (IB2b . 

For evaluation of the complicated objects like (P' |<J>[/9]E/'''(x_i_) U(y±)\ P) (see e.g. Gluon distributions) we can use a 
quasiclassical formula 



P' 



P) = 



A{x,y)B{x,y) 
d 3 Ai 

(2tt) 3 

*** f d 3 A 1 d 3 A 2 



.Ai^ 



P+ 



d 3 Xe lA1 
r f 



A| 2-2 



P 



iA2±2 



where r = x — y. 



P+ J (2tt) 3 (2tt) 3 



rf 3 A 2 

(2tt) 3 



(2 7 r) 3 (S 3 (A-A 1 -A 2 )(P + A 1 
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^'5--5 



5 + A 2 


'G- 


-0 


'» 


p\/p + A 2 


*(: 


2' 2/ 



cPX 



iAl jci 



(5-^-5-^)^(i-^-i- 



P 



(B8) 
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